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Abstract 

We study the inverse boundary value problems for the Schrodinger 
equations with Yang-Mills potentials in a bounded domain Oq C R™ 
containing finite number of smooth obstacles IX,-, 1 < j < r. We prove 
that the Dirichlet-to-Neumann opeartor on d£lo determines the gauge 
equivalence class of the Yang-Mills potentials. We also prove that the 
metric tensor can be recovered up to a diffeomorphism that is identity 
on dfto- 

1 Introduction. 

Let Qq be a smooth bounded domain in R n , diffeomorphic to a ball, n > 2, 
containing r smooth nonintersecting obstacles Qj, 1 < j < r. Consider the 
Schrodinger equation in Q = Q \ (Uj =1 Qj) with Yang-Mills potentials 

n / q \ 2 

(1.1) f -z— J m + J u + V (x)u - k 2 u = 

with the boundary conditions 

(1.2) u Uo =0, 1 < j < r, 
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(1.3) u\ dno = f(x'), 

where Aj(x),V(x),u(x) are m x m matrices, I m is the identity matrix in 
C m . Let G(Q) be the gauge group of all smooth nonsingular matrices in Q. 
Potentials A(x) = (A 1} A n ),V and A'(x) = (A[, A' n ) } V'(x) are called 
gauge equivalent if there exists g(x) G G(Q) such that 

(1.4) A\x) = g-'Ag - i9~\x)^ V = g^Vg. 

Let A be the Dirichlet-to- Neumann (D-to-N) operator on <9f2o, i.e. 



Af=(^ + i(A-v)u) 



\dn ' 



where v = (vi,...,u n ) is the unit outward normal to dQo and u{x) is the 
solution of (dUl), (JH2J), (JOJ) ). We assume that the Dirichlet problem (TTTJ) . 
()1.2|) . ()1.3|) ) has a unique solution. We shall say that the D-to-N operators 
A and A' are gauge equivalent if there exists g Q G G(Q) such that 

A' = 9 , o,an -^fl'o,9no' 

where <7o,an is the restriction of g to <9fi . We shall prove the following 
theorem: 

Theorem 1.1. Suppose that D-to-N operators A' and A corresponding to 
potentials (A', V) and (A, V) respectively are gauge equivalent for all k G 
(k Q — <5 , k Q + 5 ), where k > 0, 5 > 0. Then potentials (A', V) and (A, V) 
are gauge equivalent too. 

If we replace A', V by A® = g^A'g - ig^, V® = g^ x Vg Q then 
A = Ai where Ai is the D-to-N operator corresponding to {A^-\ V^). The 
proof of Theorem 11.11 gives that if A = A x then (A,V) and (A^\V^) are 
gauge equivalent with a gauge g G G(Q) such that g\an = I m - We shall 
denote the subgroup of G(fl) consisting of g such that g{x)\dci = Im by 
Gq(Q). In the case when Q contains no obstacles Theorem 11.11 was proven 
in [E] for n > 3 and in [E3] for n — 2. Note that the result of [E] is stronger 
since it requires that A = A^ for one value of k only. In the case n = 2 the 
proof of Theorem 11.11 is simpler than that in [E3] since it does not rely on 
the uniqueness of the inversion of the non-abelian Radon transform. 
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We shall prove Theorem II .11 in two steps. In §2 we shall prove that (A, V) 
and (A^\ V^) are locally gauge equivalent using the reduction to the inverse 
problem for the hyperbolic equations as in [B], [Bl], [KKL], [KL], [El], and 
in §3 we shall prove the global gauge equivalence using the results of §2 and 
of [E2]. Following Yang and Wu (see [WY]) one can describe the gauge 
equivalence class of A = (Ai, A n ). Fix a point x^ G dQ and consider all 
closed paths 7 in Q starting and ending at x^°K Let x = j(t), < r < r , 
be a parametric equation of 7, 7(0) = 7(70) = x^°\ Consider the Cauchy 
problem for the system 

jUr, 7 ) = ^ • A( 7 (r))c(r,7), c(0, 7 ) = I n . 

By the definition the gauge phase factor 0(7, A) is c(r ,7). Therefore A 
defines a map of the group of paths to GL(m, C). The image of this map is 
a subgroup of GL(m, C) which is called the holonomy group of A (see [Va]). 
It is easy to show (c.f. §3) that 0(7, A^) = 0(7, A^) for all closed paths 7 iff 
A^ and A^ are gauge equivalent in f2. As it was shown by Aharonov and 
Bohm [AB] the presence of distinct gauge equivalent classes of potentials can 
be detected in an experiment and this phenomenon is called the Aharonov- 
Bohm effect. In $ 4 we consider the recovery of the Riemannian metrics from 
the D-to-N operator in domains with obstacles. 



2 Inverse problem for the hyperbolic system. 

Consider two hyperbolic system: 
(2.1) 

L^u = —u^ + V(- ? ^J m + 4 p) (x))V p) + V {v \x)u^ = 0, p = 1,2, 

ot z L — ' OX a J 

j=l J 

in Q x (0, To) with zero initial conditions 

(2.2) u( p \x,0) =u ( f ) {x,0) =0 
and the Dirichlet boundary conditions 

(2.3) w W | a%x( o,To) =0, l<j<r, u w |«iox(o,2b)=/(^*). P=U- 
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Here Q = Qq \ (U r j =1 Qj) is the same as in §1, Aj p) (x), l<j<n, V®(x), 
u^(x,t), p = 1,2, are smooth m x m matrices. As in §1 introduce D-to-N 

operators A (p) J = (£ + i YTj=i A j • v i) u \an x(o,T ) , p = 1,2. 

Making the Fourier transform in t one can show that the D-to-N operator 
for (|2.1j) when T = oo determines the D-to-N operator for (|1.1JI for all k 
except a discrete set, and vice versa. 

We shall prove the following theorem: 

Theorem 2.1. Suppose = A^ 2 ^ andTo > max xg Q d(x, <9f2 ) where d(x , <9Qo) 
is the distance in from x £ Vt to dVt a . Then potentials Aj\x), 1 < j < 
n, and < j < n, V^ 2 \x) are gauge equivalent in Q, i.e. 

fT^P ao/ds a G G Q (n). 

Note that Theorem 12.11 implies Theorem 11.11 We can consider a more 
general than ()2.1j) equation when the Eucleadian metric is replaced by an 
arbitrary Riemannian metric: 



(2.4) +4 P) (a;))M (p) + y (p) (^ (p) (M) = 0, 

where lls^aOII -1 are metric tensors in Q, ,g p (x) = det A ( f\x),V^(x) 
are the same as in (j2.1|) . fi( p ) — Q \ fl , £l' p = U£ =1 % p . Let T be an open 
subset of dQo and let < T < T be small. Denote by A(0,T) the intersec- 
tion of the domain of influence of T with dQo x [0,T]. We assume that the 
domain of influence of T does not intersect Vt' p x [0, T\. 

Lemma 2.1. Suppose A^ = on A(0,T). There exist neighborhoods 
jj(p) c =1,2, TIM ndtt =T and the diffeomorphism <p : -> i/( 2 > 

stzc/i iaai y?|r = / and 1 1 #2*1 1 = V 9 ° llpi^ll- Moreover A^\ 1 < j < n, 
and <p o Aj, 1 < j ' < n, ip o are gauge equivalent in i.e. there 

exists g(x) G G (lfW), g(x) = I onT such that holds in . 

The proof of Lemma f2. 21 is the same as the proof of Lemma 2.1 in [El]. 
One should replace only the inner products of the form J u(x,t)v(x,t)dxdt 
by J Tr(uv*)dxdt where v* is the adjoint matrix to v(x, t). We do not assume 
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that matrices Aj\V {p) are self-adjoint In the latter case Lemma [2. II can be 

obtained by the BC-method (see[B], [KKL]). Extend ip -1 from {7 2 to fK 2 ) in 

such a way that p = I on 8Qq and p is a diffeomorphism of and f2( 2 ) = 

V9 _1 (fi( 2 )). Also extend gr(x) from TJ[ to fi( 2 ) so that g(x) G G (& 2 ), 9 = 1 
on dtt . Then we get that = 51 o 9? o = in U {1) . 

Lemma 2.2. Lei L 1 - 1 ) and be the operators of the form \2.4\) in QV 3 ) = 
tt \tt^, p = 1,2. Let B C n fi (2) 6e simply-connected, dB n dtt = T 
be open and connected, and D,^ \ B be smooth. Suppose = in 
B and AW = A( 2 ) on d£l x (0,T ) where are the D-to-N operators 
corresponding to L^ p \ p = 1,2. Then AW = A^ 2 ) where M p ^ are the D-to-N 
operators corresponding to in the domains (Q^ \ B) x (S,Tq — 8), 5 = 
max^gg d(x, dQo)), d(x,dflo) is the distance in B between x G B and dQ , 
are given on d(Qo \ B) x (<5, T — 5). 

Therefore Lemma 12.21 reduces the inverse problem in x (0,T ) to 
the inverse problem in a smaller domain (fiw \ B) x (8,T — 5). Combining 
Lemmas 12 . II and 12 . 21 we can prove that for any x^ G Cl^ there exist a simply- 
connected domain Bi C Vt^ l \ x^ G B\, a diffeomorphism ip of (l^ onto fK 2 ), 
^ = / on 9O , such that g G G (fiP)) such that = gopo = L« in 
To prove the global gauge equivalence and global diffeomorphism in the 
case when is not simply-connected we shall use some additional global 
quantities determined by the D-to-N operator (c.f. [E2]). 



3 Global gauge equivalence. 

In this section we shall prove Theorem 12.11 Fix arbitrary point x^ G d£lo. 
Let 7 be a path in f2 starting at x^ and ending at x^ G f2, j(t) = x(t) is 
the parametric equation of 7, < r < 77, x^°' = x(0), a^ 1 ) = x(ti). Denote 
by c^ p \t, 7), p — 1,2, the solution of the system of differential equations 

(3.1) z ^r 2 7) = ^ (r) . A (p) (x(t))c ( P ) (T)7)) 
where 

(3.2) c (p) (0,7) = ^ m , P=l,2, 0<r<n, 
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(3.3) tfr) = * 
Denote c^^'j) = (7-1,7), p= 1,2. 

Lemma 3.1. Suppose A( l > and are locally gauge equivalent. Then the 
matrix c^{x^\ 7) (cW(;rW, 7))^ depends only on the homotopy class of the 
path 7 connecting and x^ . 

Proof Let 71 and 72 be two nomotopic paths connecting x^ and x^\ 
Consider the path 7 = 71 7^ 1 that starts and ends at x^°\ It follows from 
(13. lj) that c^(r, 7)(c( 1 )(t, 7)) _1 satisfies the following system of differential 
equations: 

^(c (2 V,7)(c (1 ^ 

(3.4) -(^(^^(^(r^))- 1 )^^^)) - 7 (r). 

Let6(r, 7 ) = c( 2 )(r, 7 )(c( 1 )(r, 7 ))- 1 , b(x^, 7l ) = 6(n, 71), 72) = ^2, 72), 

where x = x^(r) are parametric equations of 7^, < r < r p , p = 1, 2. 
We have that b(x^ 1 \ , yi) = fr^ 1 -*, 72) iff b(x^°\ / yo) = I m , where x^ is the 
endpoint of path 70 = 71 72 -1 and bix^'^z) is the value at the endpoint 
of the solution of ()3.4|) along 70 with the initial value (|3.2j) . If 70 can be 
contracted to a point in f2 there exists closed paths <7i,...,er/v such that 
7o = o"i...crjv and each e>j is contained in a neighborhood t/j C f2 where 
A^ and a gauge equivalent (see Lemma 12. ip . We shall show that 
bj(r,o~j) is continuous on crj where bj(r,aj) is the solution of (|3.4j) with 7 
replaced by o~j, Cj(r) = x^'(r), < r < r,, is the parametric equation of 
°i3 = a j( T j)- The continuity on Oj means that 6(0, aj) = b(rj,cTj). 

Since and A^ are gauge equivalent in Uj there exists <7j(x) G C°°(Uj) 
such that (Ol) holds in [/,-. It follows from (JIOJ) and that 

(3.5) i^Hr^gJ^ir)) = for < r < r r 

Therefore bj{T,aj)g~ 1 {x^\r)) = C on (Tj. We have 6y(0, <Xj) = bj(Tj,o~j) = 
C gj(x^(0)). Since bj(r,cTj) is continuous on each c^-, 1 < j < N, we get 
that 6(7,70) is continuous on 70, in particular, 6(x^,0, 70) = I m . Therefore 
6(x«,7 1 ) = 6( a ;( 1 ),7 2 ). ' □ 



6 



Now we shall prove that 6(x < - 1 - ) ,7i) = b(x ( - 1 \ , j2) for any two paths con- 
necting and As in the case of Lemma f3. II it is enough to prove that 
b(r, 7) is continuous on 70 = 71 7.7 1 where b(r, 70) is the solution of (|3.4|) for 

To- _ 

We say that 7 = 71, 7^ is a broken ray in fl x [0, T ] with legs 7,-, 1 < 
j < N, if it starts at some point G dfl x [0, T ], x = x^+ru, t = 

+ r is the parametric equation of 71 for < r < T\. Then 7 makes A^ — 1 
nontangential reflections at dQ x [0,To], where Q' = U r j =1 Qj and ends at 
dfl x [0, T ]. Denote by 7 = 71... Tat the projection of 7 onto the x-plane. 
Let (r, 7) be the solution of the system 

(3.6) *|-c (p) (t,7) = #)( 7 (r)) ■ 7(t)cM(t,7), 



(3.7) c^(0, 7 )=/ m , P=l,2, 

where 7(r) is the parametric equation of broken ray, < r < r^v, 7 = ^ 
is the direction of the broken ray, c^(r, 7) is continuous on 7(1"). Note that 
-^r^ = % is constant on 7^, r^-i < r < 77 

The following lemma is the generalization of Theorem 2.1 in [E2]: 

Lemma 3.2. Let (x^ N \ t^) be the endpoint of the broken ray 7 : 7(17?) = 
(#),^)). Denote ^(^,7) = 0^(^,7), p= 1,2. Then ^'(x^j) = 
c (i)( x W ;7 ) assuming that = on <9O x (0,T ). 

Assuming that Lemma 13.21 is proven we shall complete the proof of The- 
orem 123] 

Let 7 be a broken ray starting at x^ and ending at x^ N \ x^ G dflo, x^ G 

Let c^(r, 7) be the solution of (|3.fij) . ()3.7j) . Let 01 be a path on <9fio 
connecting x^ -* and a^ 1 -* and let 02 be a path on <9f2o connecting a;^ and 
x^°\ Therefore a = 0170:2 is a closed path starting and ending at x^ . If 
Uj fl <9f2o 7^ then the gauge Oj = I m on <9f2 . Therefore Jjgy I = on 
[/j H <9f2o for any vector I G R" tangent to dflo- Then (ll.4jl implies that 

A®- 7= 7- It follows from (EH), (Q that c^(r,ai) = c( 2 )(r,oi) and 
c^^r, a 2 ) = c' 2 '(r, 02). Therefore 6(r, a) is continuous on a = 01702- We 
shall call o = 01702 an extended broken ray. 
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We shall assume for simplicity that extended broken rays generate the 
homotopy group of Q. Otherwise we can, as in the end of §2, construct a 
simply-connected domain such that fiW c ^ qq(o) d qq q ^ qq(o) n 
H 7 = and n \ is "thin", i.e. the volume of fi \ is small. Since 
is homotopic to <9f2 we get, using Lemma s 12.11 and 12.21 that potentials 
and are globally gauge equivalent in Therefore the proof of 

global gauge equivalence in Q can be reduced to the proof of the global 
gauge equivalence in Q \ Q(°\ It is clear that the extended broken rays in 
generate the fundamental group 7Ti(f2\fK°)). Note that rays without 
reflections also generated the fundamental group 7r 1 (f2\f2^°)). Then the closed 
path 7q is homotopic to where a^' are extended broken rays. Since 

b(r, ay>) is continuous on ay\ j = 1, JVi, we get that b(r, 70) is continuous 
on 7 . It follows from Lemma 13.21 that 6(r, Uj) = I m on 8Qq and hence 
6(t,7o) = I m on dil . Therefore we proved that c^^x^, ^(c^ (x^\ 7)) _1 
does not depend on the path 7 connecting x^ and x^ 1 '. Denote g(x^) = 
c^(x( 1 \ , y)(c^( y x ( - 1 \'y))~ 1 . We have that g(x) is a single-valued matrix on 
f2, g(x) = I m on dVL Q and g(x) is nonsingular since ^(x, 7) are nonsingular, 
p = 1,2. We have for arbitrary x^: 

(3.8) A^(x(r)) - 7 (r) = 7) (^(r, 7))- 1 

= ( t gg( x (r)).^T)c^+ ig (x(r))^)(c^(T, 7 ))- 1 g- 1 
d d 

= i-^g ■ i9~ l + gA (1) ■ i(r) g - 1 = i^gg' 1 + ^V 1 ) ■ 7 (r). 

Since we can choose 7(7") such that 7(71) = x^ and 7(7-) is arbitrary at 
r = n, we get that 

(3.9) A (2) (x (1) ) = i(-^g)g'\x^) + ^(xW)^ 1 )^^)^- 1 ^^), 

i.e. A( 2 ) is gauge equivalent to A^ in f2. We can change A^ to A' = 
gA^g^ + i-^gg- 1 , V = gV w g~\ Then we will have A^ = Therefore 
applying Lemmas 12.11 and 12.21 we get that = V in f2. Therefore = 
gV^g^ 1 where g(x) is the same as in ()3.9|) . □ 
It remains to prove Lemma 13.21 In the case when the broken ray 7 = 
7i..-7m does not contain caustics points the proof of Lemma l3~2l is the same 
as the proof of Theorem 2.1 in [E2]. We shall consider the case when 7 has 
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some caustics points and we shall simplify also the proof of Theorem 2.1 in 
[E2]. However in this paper we shall not use rays having caustics points. 
Consider, for simplicity, the case n = 2 and x* G Jm is the only caustics 
point on 7. We also assume that the caustics point is generic (see [V]). Note 
that if is not generic but the broken ray 7 can be approximated by a 
sequence of broken rays having only generic caustics points, then Lemma l3~21 
holds for such 7 too. This fact suggests that Lemma T3. 21 is likely true for any 
broken ray. 

Let X o(y') e C °°(R 2 ), y' = (y u y 2 ), X o(y') > 0, X o(v') = for \y'\ > 
1, Xo(y') = 1 for \y'\ < f , J R2 xl(y'W = 1. Denote 

(3.10) x (y') = \xo{^). 

We shall heavily use the notations of [E2, §2]. The difference with [E2] is that 
in this paper we consider the broken ray 7 in Q x [0,To] and its projection 
on Q will be the broken ray 7 considered in [E2]. 

Let LT be a plane in R 2 x R, (x, t) G II if x = x^ + y\0J±, t = y 2 + t^°\ 
where u>± ■ u = 0, x^ G" Cl and the plane LT does not intersect f2 x R. We 
denote by 7(2/) = 7o(?/)7i (?/')•• -1m (?/') the broken ray starting at (x^ + 
yiu±,t(°> + y 2 ) in the direction y' = (yi,y 2 ). Then the equation 

of %(y') is x = + yi u ± + t u, t = *(°) + t/2 + t , < i„ < *o(l/i)> 
where (x^ ^ + yio>j_ + £0(2/1)^3 ^ + 2/2 + ^0(2/1)) — ^1 is the point where 70 
hits dQ! x (0, To). As in [E2] we introduce "ray coordinates" (s p ,t p ) in the 
neighborhood of j p , < p < M. Denote by Dj(x(sj,tj)) the Jacobian of 
the change of coordinates x = x^\sj,tj). Let Pj be the points of reflections 
of 7(2/) at dQ f , 1 < j < M. Denote by Pj the projection of Pj on the x- 
plane. Note that the time coordinate of Pj is = t^°> + y 2 + Yli=o^r{yi) 
where t r (yi) is the distance between P r and P r -i- Note that t = + tj on 
7i, < ^ < tj(yx). 

Let Z/ p ) be the same as in (|2.1|) . p = 1,2. We construct a solution of 
L^-'u = of the form (c.f. (2.1), (2.9) in [E2], see also the earlier work [I]): 

M-l 

(3.11) u(x,t,bj)= ^ Uj(x, t, u) + U M \ + U M 2 + UM3 + U^i 

3=0 

where the principal part of Uj has a form 

(3.12) u j0 = ajofat,^*®*^-*), 
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ipj{x,u) are the same as in (2.2), (2.3), (2.4) in [E2] and 



(3.13) a j0 = (aj-. 1)0 \Dj\s) \ P . 1 Cji(x,uj), 

\ D j\* 

where Cji(x,u) is the solution of the system 

(3.14) i6j ■ Vcj-i = (A« ■ 9j) Cj i, tj^yi) < tj < tj( yi ), Cjl \ Pj = I m , 

< j < M, V = 9j is the direction of jj, 6$ = uj. 
We shall assume that 

(3.15) u (x,t,Lo) = x(y')®o 

on the plane II, i.e. when to = and x = x^ + yiLO±, t = t^°> + y%. Here «o 
is an arbitrary constant matrix. 

Let G 7m be such that is the caustics point in the x-plane. 

Note that u M1 has the same form as um-i for t < t* — Ce, where e is the 
same as in (|3.10p . solution umi is defined in a Ce-neighborhood U e of (x*, t*). 
We will not write the explicit form of um2 (see, for example, [V]) since we 
will only need an estimate 

Ckl 

(3.16) \um2\ < 



1 + kid* (x) 



where d(x) is the distance from x G Uo t£ to the caustics curve. Such estimate 
holds in the generic case (see [V]). Moreover, 

\VUM2\ < 



1 + kedi (x) 

Finally, u M3 is defined for t > £* + Ce and it has the same form as u M1 . 
The main difference is that the amplitude of %/ 3 has an extra factor e l/3 
where (3 is real. The construction and the estimate of in (|3.11|) is similar 
to [E2, Lemma 2.1] with the simplification that we consider the hyperbolic 
initial-boundary value problem with the zero initial conditions when t = 
and zero boundary conditions on dQ x (0,T ) instead of (2.9) in [E2]. Since 
we assumed that To is large enough we get that the endpoint of %/ belongs 
to dtt x (0,T ). 
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We construct a solution v(x,t,u) of L\v = similar to (|3.11|) with the 
same initial data as (|3.15|) for v with a replaced by j3 where j3 is an aritrary 
constant matrix and with the same phase function ifjj(x,u), < j < M, as 
in (ETT2J) : We have 

Af-l 

(3.17) 17= Vj(x,t,u) +v M1 +v M 2 + v M 3 + v {1 \x,t,uj) 1 

3=0 

where the principal term of v j has the following form: 

(3.18) v j0 = b j0 (x,t,u)e lk ^^- t ), 

where bjo are the same as a,jo with Cji(x,u) replaced by c*j(x, uj) where c*j 
is the solution of the system 

(3.19) i9j ■ Vcj- = (6, ■ (A®)*)^. 



Taking the adjoint of (|3.19|) we get 

(3.20) -iOrVcl^cl^eyA®). 
Denote 

(3.21) c j2 = (c;.)" 1 . 

Then fHQD implies that 

(3.22) i9j ■ Vc, 2 = {dj ■ A®)c j2 . 

We assume that satisfies zero initial conditions when t — T, x E Q, and 
zero boundary condiions on d£l x (0, T ). Substitute (|3.11J1 instead of vS 1 ' 
and (|3.17|) instead of in the Green's formula. Dividing by 2k and passing 
to the limit when k — > oo we obtain (c.f. [E2]): 

M-l „ „ T 

(3.23) °=J2 _ A(2) ) ' V ^ _ M^dt + I M , 

j=0 Jn Jo 

where Im is the integral over a neighborhood of 7a/. We make a series of 
changes of variables as in (2.43) in [E2]. 
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Note that the Jacobian D m {x^ m \sm^ m )) vanishes on the caustics set 
and therefore has a singularity there. However when we make changes 
of variables this singularity in Umi,Vmi and in um3,vm3 cancels. Note also 
that the estimate (|3.16J) implies that the integral over the neighborhood U e is 
0(yfe). Therefore taking into account that a an d flo are arbitrary matrices 
we get 

- r r 

(3.24) J2 / / X 2 (y')cj2(AM - A«) ■ e 3 c 3l dt 3 dy' + 0(y/i) = 0, 

where 7(2/) is the broken ray starting at (x^ + yiw±,t ( - ^ + y 2 ) and we use 
in (EHMJ) that = cjj (see (ET2T|) ). Note that 

(3.25) cji{AU - A®) ■ B jCjl = i8j • V{c^c jX ), 

since — cJ 2 l (A^ ■ 6 3 ) = i9j ■ Vc^ 1 . After changes of variables c 3 \ and Cj 2 in 
()3.24|) satisfy the differential equations (j3.14|) . ()3.22|1 but the initial conditions 
are different : 

(3.26) Cji \ Pj = Cj-i^p. , 1 < j < M, i = 1, 2. 

We kept the same notation for the simplicity. Taking the limit in ()3.24|) when 
e — » we get 




Since c 0i \p = I m and ()3.26|) holds we get that c M2 c mi \p m = Aru i- e - 
cmi |p m = cm2 \p m ■ Lemma E21 is proven. □ 

4 Global diffeomorphism. 

Let L^'vy) = and L^'vP 1 = be equations of the form ()2.4|) in domains 
fiW = fio \ fij,, where = U$ =1 fl ip , p = 1,2. We assume that the initial 
conditions ()2.2)l in f2*- p \ jo = 1,2 and the boundary conditions ()2.3j) with fij 
replaced by Qj P , p = 1,2, are satisfied. 
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Theorem 4.1. Suppose = A^ 2 ) on 8VLq, where are the D-to-N op- 
erators corresponding to L^ p \p = 1,2. Suppose 

T > 2 min max d v (x, dQp) 

where d p is the distance with respect to the metric tensor \\g^ Then 
there exists a diffeomorphism ip of f^ 1 ) onto fl^ such that tp = I on 8Qq 

j 1 1 jk\\ II jk\\ 

and \\g J 2 \\ = ipo \\g{ \\. 

We shall sketch the proof of Theorem 14.11 assuming for the simplicity 
that m — 1, Aj =0, 1 < J < r, p — 1,2, and T = oo. By using Lemmas 
12.11 and 12.21 we can get a simply-connected domain Vl^ C such that 
\ has a small volume. Moreover there exists a diffeomorphism (p of 

onto = p^ 1 ^ 2 ^), ip = I on 8Q such that iP^ = ip o L 2 is equal to 
T (1) in fi (0) . Note that fi (0) C H fi {2) . We also get from Lemma CO that 
A^ 1 ) = on d$}(°' \ dVlo where A^ 2 ) is the D-to-N operator corresponding 
to LS 2 > . Since f^ 1 ) is thin, there is an open subset Ti of d£l^ such that 
the endpoints of geodesies corresponding to in £1^ \ orthogonal to 
Ti, form an open subset T 2 C dVl^. Denote by Di C \ the union of 
these geodesies. It follows from the proof of Lemma [2. II (see [El]) that A^ 
on Ti uniquely determines the metric tensor ] | ^ i fc ] | 1 in the semi-geodesic 
coordinates in D\. Denote by ipi the map of D\ on Di = such that 

■01 (x) are the semi-geodesic coordinates in D\. Analogously let D 2 be the 
union of all geodesies of L^ 2 > orthogonal to T 1 and let T 2 C be the set of 
its endpoints. Denote by ip 2 (x) the semi-geodesic coordinates for iP^ and let 
D 2 = %jj 2 (D 2 ). By Lemma O^i ° L {1) =ip 2 o in Di(~)D 2 . It follows from 
Lemma O that ipj = I on I\. Note that fi (2) \ fi (0) coincide with \ fi (0) 
near Fx. 

Lemma 4.1. The following equalities hold: Di = D 2 and if) = V^Vi = 
I on T 2 . 

Proof: Since we assume that To = oo we can switch to the inverse problem 
for the equations of the form (jl.lj) . Choose parameter k e C such that the 
boundary value problem of the form (|1.1|) . ()1.2j) . (|1.3|) has a unique solution 
u p for any / E Hi(d£l^ \ OVIq) where / is the same for p = 1 and p = 2. 

Choose f nonsmooth. Denote T 2 = ipi 1 (r 2 ), r' 2 = V^Q^)- It follows from 
the unique continuation theorem that ipi o m = ip 2 o u 2 in D 2 H T^ since 
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the Cauchy data of u\ and u 2 coincide on T\. Here L^u p = 0, p = 1,2. If 
T 2 7^ r' 2 we get a contradiction since ip2 Wi is C°° outside T 2 and ^1 u 2 is 
C°° outside T 2 - Therefore Di = D 2 and ip = ip^^i is a diffeomorphism of 
Di onto D 2 . Since T 2 = T 2 we have T 2 = V'O-^)- Since ipi o m — ip 2 n 2 in 
-Di = D 2 and ui = w 2 = /(a;) on T 2 we get f(x) = f(tjj(x)) on T 2 . Since / is 
arbitrary this implies that ip — I on T 2 . □ 

Therefore ^ = / on d£>a n dQ^. Define </?W 
on £>!. We get that y?W o L< 2 ) = L^m Pu T^. 

Applying Lemma O to Q {p) \ UA) and using again Lemmas 14.11 
12.11 and 12.21 we prove Theorem 14.11 

Remark 4.1 (c.f. [El]). We shall show now that the obstacles can be 
recovered up to the diffeomorphism. Let 70 be an open subset of dtt^ 
close to the obstacle f^. Denote by Ai the union of all geodesies in f^ 1 ) 
orthogonal to 70 and ending on f^. Denote by 71 the intersection of Ai and 
Introduce semi-geodesic coordinates for in Ai. Let tpi be the change 
of variables to the semi-geodesic coordinates and let Aj = </? 1 (A 1 ). Let <p 2 
be the change of variables to the semi-geodesic coordinates for in A 2 
where A 2 is the union of all geodesies of LP orthogonal to 70 and ending 
on Q! 2 . Let 7 2 = A 2 n dQ 2 , j 2 = ^2(72), A 2 = (p 2 (A 2 ). Let = be 

a geometric optics solution in f^ 1 ) \ f^°) similar to constructed in §3 that 
starts on 70, reflects at dCl[ and leaves tt^ \ tt^ again on 7 . Let u 2 be the 
solution of = in \ having the same boundary data as U\. 

Since ipi o L^ l > = ip 2 o LP in Ai fl A 2 and since tfi o u\ and (p 2 o u 2 have the 
same Cauchy data on 70 we get by the uniqueness continuation theorem that 
(f\ou\ = ip 2 o u 2 in Ai fl A 2 . If 71 7^ 7 2 then we can find u\ such that tpi o m 
and ip 2 o u 2 will have different point of reflection and this will contradict that 
ifi o u\ = ip 2 o u 2 in A x fl A 2 . Since 71 = 7 2 we get that <£>(tl) = 7 2 C dQ! 2 
and ip(Ai) = A 2 where if = (p^fi- □ 

Remark 4.2 Note that Lemma 14.11 allows to consider the inverse prob- 
lems in multi-connected domains fl with the D-to-N operator given on a not 
connected part r of dfl. 
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